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1 Abstract 

We draw attention to a problem with the alpha-Omega dynamo when it is 
applied to the origin of the galactic magnetic field under the assumption 
of perfect flux freezing. The standard theory involves the expulsion of un¬ 
desirable flux and, because of flux freezing, the mass anchored on this flux 
also must be expelled. The strong galactic gravitational field makes this im¬ 
possible on energetic grounds. It is shown that if only short pieces of the 
undesirable field lines are expelled, then mass can flow down along these 
field lines without requiring much energy. This expulsion of only short lines 
of force can be accomplished by a spike instability associated with gigantic 
astrophysical superbubbles. The physics of this instability is discussed and 
the results enable an estimate to be made of the number of spikes in the 
galaxy. It appears that there are probably enough spikes to cut all the unde¬ 
sirable lines into pieces as short as a couple of kiloparsecs during a dynamo 
time of a billion years. These cut pieces then may be randomly rotated in a 
dynamo time by alpha-Omega diffusion and there is enough rotation to get 
rid of the undesirable flux without expelling the fields themselves. The spike 
process seems strong enough to allows the alpha-Omega dynamo to create 
the galactic field without any trouble from the boundary condition problem. 


2 Introduction 

Our galaxy is believed to have hnite magnetic flux in the form of a toroidal 
field, because this field does not reverse across the galactic midplane. The 
origin of this field is paradoxical. This is because, on the galactic scale, flux 
freezing is almost infinitely strong so that flux through any moving plasma 
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region cannot change. In fact, the change is only hnite over a Hubble time if 
the scale of variation is smaller than an astronomical unit. Thus, the obvious 
question arises: How can one start with a weak held of cosmological origin 
and increase it to its present value? 

The answer is: that the hux in the galactic disc itself need not be constant, 
but only the hux in the larger region consisting of the disc plus the galactic 
halo. The standard alpha-Omega theory (Steenbeck et. al [1]) supposes one 
starts with a very weak held whose origin is cosmological, and is amplihed 
by compression during galactic formation. (Such an initial held might have 
a strength of Bq = 10“^^ gauss, be toroidal in the positive direction, and 
uniformly hll the disc with hux $.) The standard alpha-Omega dynamo 
folds this held back and forth by the alpha ehect. The result is a hux of 2$ 
in the positive direction about the midplane and a negative hux of —<h near 
the edges of the disc. Then the negative hux is supposed to be turbulently 
dihused out of the disc and into the halo leaving the disc with double its 
original hux. This doubling takes place during a ’dynamo time’ of less than 
a billion years. It repeats over and over perhaps twenty times in the life of 
the galactic disc amplifying the held strength (and hux in the disc) by over 
a million times. Such a process appears to provide a reasonable origin of our 
present galactic held, and also satishes the magnetic hux freezing condition, 
(Parker [2], Ruzmaikin, et. al. [3]). 

However, when one looks more closely at the dihusion of hux from the 
disc to the halo, a problem arises because the galaxy has a very strong 
gravitational held. In the motion of a hux tube into the halo, the mass 
anchored onto it by hux freezing must also be lifted into the halo. 

During this motion an energy, equivalent to an escape velocity of four 
hundred kilometers a second, must be supplied to the mass. This energy is 
extremely large compared to energies in interstellar turbulence, which have 
typical velocities of ten kilometers a second. Where can such energies come 
from? 

One suggestion is supernovae. But velocities of supernovae remnants 
are reduced to a few tens of kilometers a second by snow plowing all the 
surrounding interstellar mass. They have only a small chance of breaking 
out of the galactic disc. 

An alternative suggestion is the recently recognized phenomena of super¬ 
bubbles, which are driven by multiple supernova. Some of these actually 
break out of the galactic disc. However, by the time they leave the disc, 
the mass that they have snow plowed has slowed them down to a velocity of 
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fifteen to twenty kilometers a second. So this can not be a direct source of 
the large energies needed for diffusion of the tubes into the halo. 

One is forced to consider a different concept. This concept is based on the 
realization that the flux freezing process applies only to motions perpendicu¬ 
lar to the line of force. It in no way constrains the parallel motion along the 
held. Thus, if a short length Is of a line of force is lifted out of the disc and 
into the halo, the mass anchored on it can be greatly reduced from its initial 
value by sliding down along a length ^v to the rest of the line remaining 
in the disc. Because of the small mass on the ^s little energy is required. 
Thus, on energy grounds alone, such a process is possible. But removing a 
short piece of a negative line into the halo does not itself provide a reduction 
in the negative toroidal hux in the disc. A further idea is needed (Kulsrud 

[4]). 

Consider what happens to the rest of the negative held line on either 
side of £ 5 ? There is a gap in it caused by the removal of £ 5 . Of course, 
this gap is kept closed by two vertical pieces of the held connected to Ig- 
But the held strength of these pieces, t'y, becomes very weak by horizontal 
expansion. The region into which these expand is much larger than the 
thickness of the disc by a factor of hfty, the aspect ratio of the disc, so that 
the held strength of these ^v become far too weak to ahect the rest of the 
line, Ir. This part of the line acts as though its two pieces are ehectively 
decoupled or cut. Therefore, the ends of Er are free to move and it appears 
that this part of the held line has free ends. 

Suppose a given line of force is ’cut’ into a number of pieces by removal 
of a number of short is pieces. Then, when these Ers are acted on by (3, 
the turbulent dihusion of the alpha-Omega dynamo, they will be rotated 
randomly in direction and no longer preserve their negative toroidal hux. 
The rotations transfer their negative hux to the is pieces in the halo. Due to 
the weakness of the iy pieces this transfer through the weak iy pieces of hux 
will have no dynamic ehect on the disc pieces iR. But this transfer ehectively 
removes negative hux into the halo as required by the alpha-Omega dynamo. 
Only a small amount of energy is required for these processes. See hgure 1. 


It remains to discover a mechanism to take the small is pieces and propel 
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Figure 1: (a) A short piece of a field line is is lifted into the halo. It is 
connected to the pieces in in the disc by weak pieces iy- (b) There is a gap 
in in that allows the ends of the held lines to rotate freely under diffusion so 
the points a and b rotate to new positions. 
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them far into the halo. In paper I (Kulsrud [5]) it is proposed that super¬ 
bubbles present a way to do this. 

Each superbubble arises from the formation of a normal galactic star clus¬ 
ter. In the star cluster the stars are formed in a short time and the more 
massive O and B stars in the cluster continually evolve and explode into su¬ 
pernovae. The many supernova explosions overlap and their energy is trans¬ 
formed into a high pressure low density plasma. Its pressure forces all the 
surrounding interstellar matter into a gigantic rapidly expanding shell which 
comprises the superbubble phenomena. The shell gradually slows down as all 
the surrounding mass accumulates onto it. The shells of the more luminous 
superbubbles manage to escape out of the disc. Once they leave the disc, 
their expansion accelerates because the pressure of the core remains large 
(the supernova are still exploding), and there is no more mass to snowplow. 

This situation where low density plasma is accelerating the high density 
shell is unstable to the Raleigh-Taylor instability, (MacLow and McCray [6]) 
and the shell breaks up into many fragments. The fragments may still be 
connected and the pressure still contained by the fragmented shells. This 
leads to a secondary instability the ’spike instability’. 

This instability is described in paper I. The spike is a narrow elongated 
eruption on a fragment. The matter at the top of the spike is accelerated 
against gravity, and because of this gravity the mass falls down along the 
sides of the spike. (See hgure 2). As the mass at the top slides down the 
top becomes lighter and, as long as the pressure is confined, it accelerates 
faster. This increases the steepness of the sides of the spike and the down 
flow increases, making the top still lighter and accelerate faster. This process 
I call the ’spike instability’. As the spike rises towards the halo it increases 
its speed until hnally the velocity exceeds the escape velocity and the top 
mass reaches the halo ballistically. 


The shell mass contains the interstellar field and as the spike rises from 
the fragment, it carries the is pieces of the same lines of force with it. Most 
of the lines slide down with the mass but some of those at the very top of 
the spike reach the halo as desired for the model of the dynamo. 

For each spike the i^ parts of these lines of force are cut. When enough 
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Figure 2: As the spike grows, gravity g leads to a downflow of the mass from 
the top of the spike to the fragment, making the top lighter. The superbubble 
pressure then lifts the top faster 
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lines are cut into short enough pieces the negative flux can be diffused away 
by random turbulent rotation of these Ir pieces. If this occurs during a 
short enough time compared to dynamo times (say a billion years), then the 
alpha-Omega dynamo should work and the galactic held amplified from a 
weak held to its present value. 

3 The Spike Instability 

In this section we develop the theory of the spike instability in order to 
determine what fraction of the magnetic lines of force initially in the spike 
are expelled to the halo. 

Then we estimate the number of spikes and the average length of the 
cut lines, and Anally, we calculate the amount of angular randomization of 
the lines that occurs during a dynamo time. It from this calculation it 
appears that the process of flux removal works during a dynamo time of a 
billion years. In the final section, we summarize and draw the conclusions. 

As a model for the spike we approximate the surface of the superbub¬ 
ble fragments as infinitely thin plane sheets with surface mass ao supported 
against gravity by a constant underlying pressure. Its surface density and 
actual thickness D are derived in the paper I. The gravitational acceleration 
g is the sum of the true gravity due to the stars in the disc and the effec¬ 
tive gravity from the acceleration of the superbubble shell. However, in our 
examples the astrophysical gravity is the more important. 

We work in cylindrical coordinates with z vertical to the sheet and r and 9 
in the initial plane of the sheet. We write the equations of the perturbed 
motion in a Lagrangian form. We assume cylindrical symmetry and assume 
all quantities depend only on time t and the initial radius of the perturbation 
tq. Let ^r(t,ro) and ^z(t,ro) be the components of the displacement. 

Let the coordinates of the four corners of an undisplaced small square Aq 
at time t = 0 be 


PI = 

[ro, 0, 0] 

P2 = 

[ro -l- 5, 0, 0], 

P3 = 

[’^o, e, 0] 

P4 = 

[ro -h S, e, 0] 


( 1 ) 


Then after the displacement, the coordinates of the four corners of the 
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displaced square A are 


PI = [ro+ Cr(?^o),0,C2(^o)] : 

P2 = [ro + <5 + e.(ro) + 5%^,0,e.(ro) + <5^|^; 

d^zjro) . 

drn 


( 2 ) 


dro 

P3 = [ro+ Cr(?^o),e,G(’^o)] 

P4 = [ro + 5 +'Cr(^o) + <^—'Cz(^o) + <^ 


For notational convenience, let us abbreviate the components of the dis¬ 
placement as 


^r{ro) = R (3) 

e.(ro) = ^ 


After the displacement, the components of the vector A are 

Ar = -e5rZ' = -e5{ro + R)Z' (4) 

= e6r{l +R’) = e6{ro +R){1 +R') 


where prime denotes a derivative with respect to tq. (Where convenient we 
drop the arguments t, and tq). The vector A is perpendicular to the displaced 
square and its magnitude is equal to its area. 

In the equation of motion of the surface of the spike the mass of the 
square is (To|Ao|, the pressure force on it is pA, and the gravitational force 
is —gao\Ao\z. From the initial equilibrium we have gaQ = p 
Thus, the equations of motion of the displacement are 


d^R 

dR 

d^Z 

dR 


-9(1 + -)Z' 

g{l H-)(1 + R!) — g 

ro 


(5) 


where the prime denotes the derivative with respect to tq. 

For the linear solution of these equations assume that R and Z are pro¬ 
portional to exp 'jt. Keeping the linear terms we have 


= -gZ' 

= =g(- + R') 

ro 

8 


( 6 ) 



Defining 


(7) 


k = 


7 


we find 


kg 

Z” + — + k^Z 
ro 

i?" + :^ + _ 4 = 0 

ro ri 


= 0 


so the linear solutions regular at ro = 0 are 
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Z = ae^KJo{kro) (9) 

R = ae^*Ji{kro) 


The growth rate is 7 = \/%. The shortest wave length consistent with an 
equilibrium sheet mass is given by kD = 1 where D is the actual thickness 
of the layer which we approximated above as a sheet. An estimate of this 
thickness, which is given in paper I, is about a parsec. The value of g in the 
neighborhood of the disc is 10“^ cm/ sec^. These values give the linear time 
of growth, fc = 1/7 = 7.4 x 10® years, a time short compared to the life time 
of the superbubble (usually 20-50 million years). 

To get the long time behavior of the spike, we approximate Equation ([5]) 
by its nonlinear terms. 


d^R 

w 

d^z 

dR 

These equations have a solution 
Setting 


= -—RZ' 

To 

= ^RR! 
rQ 


( 10 ) 


with R and Z proportional to l/(fo — ^)^- 


R 

Z 

We have from the hrst equation. 


V 

{to - ty 

V 

(V^ 


67 



( 11 ) 


( 12 ) 
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( 13 ) 


Canceling rj from this equation, and solving for z/ gives 


The equation for rj is 


u = C-3-^ 
9 


6 z/ = — 7777 ' 
^^0 


which on substituting the solution for u gives 

, 12 c 36 3 

(V ) = ^ro - ^ro 


2 6C 2 94 6C 2 

V = —ro - ^ —Tq (16) 

The constants C and to are determined by matching both terms of the 
solution for Z to its linear solution at t = 0. 

= „(i _ erin) (17) 


which gives 


C 

72 “ ® 

bn 


Solving these two equations for to and C we have 


to — 


^ - “‘-FJ 


12 (1 - k'^rl/A) 


k‘^9 {to - ty 

Treating R the same way we get 


R== 


kg {to - tf 
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Now let us determine the time, , at which the upward velocity of the 
peak of the spike reaches the escape velocity from the of the disc, Ve, (Heiles 
[7]). That is, from dZ/dt = Ug, (Equation fl20|) h 


we get 


dZ _ 24 

dt k‘^g{At)l 


{At), = t-t. 



This occurs at the height 


hg 


1.44 



1/3 


( 22 ) 


(23) 


(24) 


With the above numbers h, = 49 parsecs. 


4 The Number of Lines Cut by One Spike 


To determine the number of lines expelled we start with the number of lines 
initially embedded in the spike. This is equal to the number of lines embed¬ 
ded in the superbubble shell in a strip with the width of one spike, ^ X = 1/k 
Since the shell consists of all the interstellar mass swept up by the superbub¬ 
ble, the strip contains BXH lines, where B is mean held in the disc. At the 
escape time tg, the region ro < A has horizontally expanded a distance 


R{te,X) ^ 1 

ro kg {At)l 


72V2 /kvlV^^ 
( 24 ) 2/3 ^ g J 


37.59pc. 


(25) 


This expansion reduces the held strength (i.e. the density of the lines) at the 
top of the spike, by a factor of 


^0 

R(t„ A) 


.0266 


(26) 


However, by the time the spike has reached a height of h,, the solution 
of the spike has spread out a distance of 37 parsecs, and the behavior of the 
spike a distance r S> A from the cylindrical axis is uncertain. Therefore, we 
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take for the number of expelled lines, <he, the number of lines at the top of 
one spike, within tq < tX. (Those that are certainly expelled.) 

= 2 X .0266 X BHtX (27) 


The other lines may fall back down, or possibly be expelled. 

Among the these lines there are both positive and negative lines. The 
proportions depend on how high above the galactic midplane the superbubble 
started. Also, the sheet of plasma, assumed in our derivation to be inhnitely 
thin, actually has a hnite thickness. In the radial outflow flow, the mass 
in the sheet nearest the superbubble pressure will probably be squeezed out 
faster than that farther away. Without analyzing this further we simply 
assume that half the number of lines at the top of the spike are negative and 
ignore the positive lines. Thus, we take the number of negative lines cut per 
spike to be 

= .0266BHXt (28) 

Denote the number of spikes per fragment as fs, and the number of 
fragments per superbubble as Nf. The size of the spike is A ~ D where D is 
the thickness of the superbubble shell. 

Consider a superbubble of luminosity L = 1.37 x 10 ^^ ergs per second. 
This is the luminosity of a superbubble just strong enough to break out of 
the disc, and fragment. For such a superbubble we have, from paper I, that 


Nf 

D 


78 

-'300 


H 


-T 


300 


(29) 


where the temperature, T = 300 x T 300 degrees Kelvin, is the temperature 
in the superbubble shell. T is generally assumed to be between a hundred 
and a thousand degrees Kelvin. Then our estimate for the number of cuts of 
negative lines by a single superbubble ^sb is 


^sb = ^^eNffs = .0266BDHNffsT 
Substituting the numbers from Equation fl29|) 


^sb — 


.0266 


78 BHy-sr 
222 T300 


0.00935 


T300 


(30) 


(31) 
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Now what is the number of lines cut by a superbubble with luminosity 
L > LI Again from paper I we hnd that Nf {L/Lf/^ and D (L/L)‘/3 
SO the number of negative lines cut by a luminosity-L superbubble is 


= ^sb 



(32) 


Consider an annular region in the galaxy with mean radius Rs = 8.5 kpc 
(the solar galactic radius) and radial thickness AR. The rate of birth of 
superbubbles with luminosity in dL is (Ferriere [8]). 

da = 0.86 X 10“^ kpc“^ yr“^ (33) 

and the rate of birth in the entire annulus is 


2TiRsARda 


(34) 


so that the rate of line cutting by all the superbubbles with luminosity L in 
the annulus is 

^Lda2TiRsAR (35) 

Now, the number of magnetic lines of force in the annulus, assuming the 
lines are toroidal, is 2BHAR. But only a quarter of these lines are negative, 
so their number, is 

V’B = .5BHAR (36) 

The rate of cutting of any given line pc by all the superbubbles with 
luminosity greater than L is the rate of cutting by all these superbubbles 
(integrated over their luminosity) divided by the number of lines, ips- 


Pc 


2'kRsAR ^Lda^L 
ABHAR 

An X .00985 x (.86 x 10-'^)RsH 
An X .00935 x (.86 x 10"^) x RsH 


dL - 2 . 3 - 1/3 

Li T [l 

fsr 


^30o(1-3 + 1/3) 


.17 X 10 


fsr 


-7 


fsr 


T3oo(1.3 + 1/3) 


yr 


-1 


10.5—— cuts per billions years 

Tqi 


^300 


(37) 
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Thus, in a billion years, each line of length 2tiRs = 53 kiloparsecs will be 
cut into 10.5/5r/T3oo pieces of lengths I = STTsoo/Zst kiloparsecs. 

These lines will be rotated through an angle by the same turbulent 
diffusion (3 as included in the normal alpha-Omega dynamo. The rate of 
spread in angle is 


{Aef _ f3t 
2 “ ¥ 


(38) 


Take /9 = /926 x 10^® cm^/ sec. I32g ~ 0.5 (Ferriere [8]). Then 


AO 


0.52 






kpc 

fsT 


T. 


300 


(39) 


where fg is the time in billions of years and ^kpc is the length of the cut pieces 
of lines in kiloparsecs. 

This, estimate indicates that the boundary condition on the alpha-Omega 
dynamo may still be a problem. The uncertainties of the estimate are en¬ 
capsulated in the parameters /s',r, and Taoo- If the temperature T in the 
superbubble shell were as low as 100 degrees, and the parameter r was as 
large as 3, then the random rotation 9 would be of 0.9. This value for 9 
is large enough to satisfy the alpha-Omega galactic dynamo with a dynamo 
time of a billion years. This choice of the parameters is quite reasonable. 

It should be noted that, in this model for flux expulsion, the galactic 
magnetic held has been assumed weak enough that any magnetic tension 
is too weak to interfere with the spike instability. This limit on the held 
is reasonable since the main problem in the origin of our galactic magnetic 
held is how to amplify it when it is extremely weak. When the galactic 
held reaches its present value, presumably this back reaction of the held will 
saturate the hux expulsion process and end any further growth of the held. 


5 Conclusion 

The purpose of this paper is to form a more precise estimate of the hux 
expulsion by spikes associated with superbubbles, than was given in paper 
I. The nature and physics of the spikes certainly introduces considerable 
uncertainty into the conclusions as to whether superbubbles can produce the 
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required destruction of negative flux needed for the dynamo to act on a time 
scale of a billion years. In this note we present a conservative estimate of the 
superbubble process that shows that one should, take the spikes seriously 
in dynamo theory. A multi-dimensional simulation of the spikes is clearly 
necessary to properly evaluate whether this process is a reasonable way to 
make the alpha-Omega dynamo viable for our galaxy. However, it should be 
noted that there are serious difficulties with the dynamo itself, and at the 
moment there does not seem to be any other reasonable alternatives. 
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